Thin front propagation in steady and unsteady cellular flows 
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Front propagation in two dimensional steady and unsteady cellular flows is investigated in the 
limit of very fast reaction and sharp front, i.e., in the geometrical optics limit. In the steady case, 
by means of a simplified model, we provide an analytical approximation for the front speed, Uf, as 
a function of the stirring intensity, U, in good agreement with the numerical results and, for large 
U, the behavior Vf ~ U/ log (U) is predicted. The large scale of the velocity field mainly rules the 
front speed behavior even in the presence of smaller scales. In the unsteady (time-periodic) case, 
the front speed displays a phase-locking on the flow frequency and, albeit the Lagrangian dynamics 
is chaotic, chaos in front dynamics only survives for a transient. Asymptotically the front evolves 
periodically and chaos manifests only in the spatially wrinkled structure of the front. 
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I. INTRODUCTION 

Front propagation in fluid flows is relevant to many 
fields of sciences and technology ranging from marine 
ecology to chemistry HQ] and combustion technol- 
ogy [j5| . A complete description of the problem would re- 
quire to consider the coupled evolution of reactants and 
velocity field, including the modification of the advecting 
field induced by the reaction. In general, this is a very 
difficult task ||. Here we consider a simplified, but still 
physically significant, problem by neglecting the influ- 
ence of the reactants on the velocity field. This amounts 
to consider the reaction as a constant-density process. 
Aqueous auto- catalytic reactions, and gaseous combus- 
tion with a large flow intensity but sufficiently low values 
of gas expansion across the flame are important examples 
of chemical-physical systems for which this approxima- 
tions is appropriate 0. 

In the simplest model a scalar field, 0(r, t), which repre- 
sents the fractional concentrations of the reaction's prod- 
ucts (i.e., 9 = 1 indicates inert material, 8 = fresh one 
and < 8 < 1 means that fresh material coexists with 
products), evolves according to the advection-reaction- 
diffusion equation 



d t 8 + u ■ V8 = D A8 + -f(8) , 

r 



(1) 



where Dq is the molecular diffusivity, and u is a given 
incompressible (V • u = 0) velocity field. f(8) is the 
production term and r its time scale. 

Two limiting cases of Eq. (Q) have been previously 
studied in detail: f(8) = and u = 0. In the former, the 
equation for a passive scalar is recovered (for a review 
see Ref. flCj). The latter corresponds to the reaction- 
diffusion equation, which has gathered much attention 
since the seminal works of Fisher and Kolmogorov- 
Petrovsky-Piskunov (FKPP) 0^2) (see also Ref. j| and 
references therein). 




Eq. (y) can be studied for different geometries and 
boundary conditions. For instance, one can consider an 
infinite strip in the horizontal direction with a reservoir 
of fresh material on the right and inert products on the 
left, and periodic boundary conditions along the trans- 
verse direction. With this geometry a front of inert ma- 
terial (stable phase) propagates from left to right. If 
the medium is at rest with the FKPP production term, 
f(8) = 8(1 — 8), the front propagates with an asymptotic 
speed and thickness given by pjll],|12| 



(2) 



where c is a constant depending on the definition adopted 
for £. This result is valid whenever f(8) is a convex func- 
tion (/" < 0) with /(0) = /(l) = and /'(0) = 1. For 
non-convex f(6) only upper and lower bounds for the 
front speed can be provided jj| . 

A more interesting physical situation is when the veloc- 
ity field is non zero. In this case, generally, the front prop- 
agates with an average limiting speed, V{, enhanced with 
respect to the fluid at rest (vf > Vq)- For very slow re- 
action, by means of homogenization techniques jlCj], one 
can show that the front speed behaves as in Eq. (pp with 
D replaced by a renormalized diffusion coefficient, D c h, 
the so-called eddy diffusivity (see Ref. [|l0| for an exhaus- 
tive review on the determination of the eddy-diffusivity). 
In realistic systems, the reaction time scale is of the same 
order or (more frequently) faster than the velocity time 
scale (fast reaction), so that a simple renormalization of 
Dq is not sufficient to encompass the dynamical proper- 
ties of the system . In some cases the front speed can 
be still obtained by Eq. @ renormalizing not only the 
diffusion constant, Dq — > D e ff, but also the reaction time 
t — > T c ff |Q. However, a general method to compute Vf 
for a generic velocity field does not exist. 

Here we consider the limit of very fast reaction and 
very thin front, i.e., the so-called geometrical optics 
regime [jj] . Formally, this corresponds to the limit r — > 



1 



and Dn — » maintaining the ratio Dq/t constant p5[ : 
from (g) this means that i>o is finite and £ — > 0. In this 
regime the front is identified as a surface (a line in 2d), 
and the effect of the velocity field is to wrinkle the front 
increasing its area (length in 2d) and thereby its speed || . 

As to the velocity field we consider steady and un- 
steady cellular flows (i.e., with closed streamlines) in two- 
dimensions p6|-pj| . Since coherent vortical structures are 
typically present in real hydrodynamical systems, cellular 
flows offer an idealized (but non-trivial) model to study 
the effects of this kind of structures on front propagation. 
Real flows, e.g., turbulent flows, are usually character- 
ized by a very complex temporal dynamics and spatial 
development of scales. In this respect a steady cellular 
flow is too simple. Therefore, we also consider either the 
presence of small scale spatial structures in the velocity 
field, or the effects of time dependence, which induces 
a complex temporal behavior for particle trajectories - 
Lagrangian chaos . 

We found that the front speed is mainly determined by 
the large scales velocity characteristics, i.e., the effect of 
small scales is just to renormalize the largest scale veloc- 
ity intensity. Concerning time periodic cellular flows, we 
found that the front speed is not significantly modified 
with respect to the steady case, although rather subtle 
effects appear. Namely the front speed locks on the flow 
frequency: a phenomenon which goes under the name of 
mode- locking (or frequency-locking) |2^,|2^], and which 
has been already found in some models of front propaga- 
tion p4[ . Lagrangian chaos is suppressed by the reaction 
and only survives for a transient. However, asymptoti- 
cally the spatial wrinkling ( "complexity" ) of the front is 
enhanced with respect to the steady case. We introduced 
a suitable observable to quantify the front complexity. 

The paper is organized as follows. In Section [n] we dis- 
cuss the geometrical optics limit. Numerical results for 
steady cellular flows with one and more scales are pre- 
sented in Section III, where we propose a simple model 
which well reproduces the numerical results. The effects 
of Lagrangian chaos and velocity phase-locking in time 
dependent cellular flows are discussed in Section Fi- 
nal remarks are reported in Section The Appendices 
are devoted to the numerical methods here employed and 
to a more detailed treatment of the frequency-locking 
phenomenon. 



II. THE GEOMETRICAL OPTICS LIMIT 

From a physical point of view the geometrical optics 
limit (in combustion jargon, the flamelet regime) cor- 
responds to situations in which the reaction time scale 
and reaction zone thickness are much faster and much 
smaller than the time and length scales of the flow, re- 
spectively (e.g., in turbulent flows this means that the 
front thickness is smaller than the Kolmogorov length 
scaled, §• 



Being the front sharp, its dynamics can be described 
in terms of the evolution of the surface (line in 2d) which 
divides the inert material (9 — 1) from the fresh one 
(9 = 0). In this limit the problem can be formulated in 
terms of the evolution of a scalar field, G(r, i), where the 
iso-line (in 2D) G(r, t) — represents the front: G > is 
the inert material and G < is the fresh one. G evolves 
according to the so-called G-equation flslJIf 



dG 
~dt 



u- VG = v Q \VG\ 



(3) 



The analytical treatment of this problem is not trivial, 
and even in relatively simple cases (e.g., shear flows) nu- 
merical analysis is needed. 

Recently Majda and collaborators |^7| pointed out 
that there are situations in which the G-equation fails 
in reproducing the front speed of the original reaction- 
advection-diffusion model. Indeed, in some systems the 
exact treatment of Eq. (]l|) in the limit r — > 0, Dq — > 
with Dq/t = const does not lead to the same results 
of the G-equation. However, for the application we are 
interested in, the study of the G-equation is physically 
relevant |l8[ . 

In the absence of stirring (u = 0) the front evolves ac- 
cording to the Huygens principle, i.e., a point x belonging 
to the front surface moves with a velocity v(x) = voh(x), 
where n(x) is the perpendicular direction to the front 
surface in x. With open boundary conditions, at large 
times the front surface is asymptotically close to a sphere 
(circle in 2c?). However, the preasymptotic behavior is 
mathematically non trivial |29| and interesting in some 
technological problems. 

In the presence of stirring (u ^= 0) the problem is much 
more difficult. The first attempt to determine the front 
speed in such a regime dates back to the 40's with the 
work of Damkoler Q who suggested that, if the veloc- 
ity field does not change the local (bare) front speed, i>o, 
then the effective front speed is proportional to the to- 
tal front area divided by the cross-section flow area. In 
two-dimensional geometry this means that 



Vf/v = Lf/L 



(4) 



where L is the transverse length, Vf and Lf are the aver- 
age front speed and length respectively. 

The average velocity and length are V[ = (v(t)) and 
Lf = (C(t)), where (■) indicates a time average. The in- 
stantaneous front velocity, v(t), and length, C(t), can be 
defined as follows. The velocity is given by flq] 



v(t) = d t {^J^v 



dx 9(x, y; t) 



(5) 



In order to define the instantaneous front length, C(t), 
we introduce the variable a t (x,y;t) which assumes the 
value if 9 is constant inside a circle of radius e centered 
in (x,y), otherwise cr t (x,y;t) = 1 (i.e., a e (x,y;t) = 1 
only if the e-ball centered in (x,y) contains a portion of 
the front). The front length is then defined by 
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C{t) 



lim — 



Ax 



(6) 



Vf 



lim 

t — >oo 



X M (t) 



(9) 



III. STATIONARY CELLULAR FLOW 

We consider the following two-dimensional cellular 
flow, originally introduced in Ref. f30|| to mimic roll pat- 
terns in Rayleigh-Bernard convection, 



Usm(^-x) cos(^y) 



u x (x,y) 

u v {x, y) = -U cos (^-x) sin (^y) 



(7) 



where [/ is the flow intensity, L the roll size (in the follow- 
ing for simplicity, L — 2n) . Periodic boundary conditions 
in the transverse directions are assumed. By fixing 9 = 1 
for x — ► — oo and 9 = for x — > oo the front propa- 
gates from left to right. It is natural to expect that the 
front speed can be expressed according to the following 
relation RSF 



Vf = 



(8) 



Our aim is to investigate the dependence of vt on the flow 
intensity U. 




FIG. 1. Front speed as a function of time, measured in the 
standard way v(t) (B (□) and as xnj(t)/t (|^)(o). The straight 
line is the average front speed Vf. The system parameters are 
U = 4, vo = 1 and L = 2ir. 

Due to the spatial periodicity of the flow, after an ini- 
tial transient, the front propagates periodically in time 
(with period T) . In Fig. [j] a typical time series of the in- 
stantaneous velocity v(t) is reported: peaks occur when 
the front length is maximal. The front speed can be de- 
fined as the average Vf = (u(£))t over a period. Since the 
interface is sharp, we can easily track the farther edge of 
the interface (xm(£), 2/m(£))> i- e -j the rightmost point (in 
the x-direction) for which 6{xm, Vm', t) — 1. Thus we can 
define a velocity 



which, since the front dynamics is periodic, is equivalent 
to the standard definition (^) (see Fig. |l|). 

In Fig. H we show some snapshots of the front at dif- 
ferent times. It is possible to see cusps and pockets of 
unburnt material (white) left behind the front edge, as 
firstly noticed by Ashurst and Sivanshinsky |H| . At high 
field intensity a trail of pockets is formed. 

We now revert to the behavior of Vf as a function of U. 
As far as we know, apart from very simple shear flows (for 
which ip(U) = U + 1) |32|, there are no general methods 
to compute ip(<U) from first principles. For the cellular 
flow under investigation, by mapping the front dynamics 
onto a one dimensional problem, it is possible to obtain 
an approximate expression for ip in good agreement with 
numerical data. 



4 




r 






FIG. 2. Snapshot of the front shape with time step T/8 
(from (a) to (d)), where T is the period of the front dynam- 
ics, for vo = 0.5, U = 4.0 and L = 2tt. 

The strategy is to devise an equation for the edge point 
evolution (shown in Fig. ||). After the transient, in the 
cell [0,2n] x [0, 7r], the point (xm(£), J/m(£)) approxima- 
tively moves in the right direction along the separatrices, 
so that yM(t) practically assumes values close to or 
7r. Therefore, one can reduce the edge dynamics to the 
following ld-problem 



d.r 



A I 



dt 



v Q + Uf3\ sin(x M )| • 



(10) 



The second term of the r.h.s. mimics the horizontal 
component of the velocity field, where /3 takes into ac- 
count the "average" effect of the dependence on the ver- 
tical coordinate, y. The periodicity of the r.h.s. of 
Eq. (0) implies that there exists a time Tm such that 
XM(t + Tm) = XM{t) + 7T. Then, by solving ( |l0| ) in the 
interval xm € [0, n] we computed Tm and hence the front 
speed, vt = n/T M , 



2 In (U(3+y/(Upy 



(11) 
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Notice that (|Tl]) is valid only for Uf3 > 1, which is the 
regime we are interested in. Now the problem is to esti- 
mate (3. 

As stated previously the front evolution is periodic, so 
that yM{t) is a periodic function of time with period T y 
commensurate to Tm (see Fig. 3); the front period T 
should also be commensurable to Tm and T y . Therefore, 
for a specific value of U and vq we numerically identified 
the period T y and computed f3 as 



cellular flow here studied, Aldredge |25],^6| compared his 
numerical results with (O) for a = 2 and d = 1/ \[2 find- 
ing a fairly good agreement (our data and Eq. (O) are 
shown in Fig. ||). Actually the asymptotic behavior of 
Eq. ( fl4| ) is Vf ~ U I '(h\U) l / a which, by comparing with 
ours results, would suggest a = 1. In the following we in- 
vestigate how small spatial scales modify the dependence 
of Vf on the flow intensity. 



1 



T, 



13 =— |cOB(tf M (i))|df, 



(12) 



y Jo 



obtaining (3 ss 0.875. Using this value in ( |TT| ) one has a 
remarkable agreement with the measured Vf, see Fig. |] 
(the agreement is between 6% and 10% for all the range 
of investigated values of U and v ). Notice that, by def- 
inition, (3 < 1. Therefore, with [3 = 1 in Eq. ( |ll|) one 
obtains an upper bound for the front speed. Moreover, 
in Ref. B3j, a rigorous lower bound has been provided: 



vt>U I Hl + U/vo). 



(13) 



It is worth remarking that, for large U, both the upper 
and lower bounds give Vf ~ U / In U, which identify the 
asymptotic behaviour of the front speed. 




FIG. 3. Time evolution of the edge point: xmW) and j/m(£). 
The simulation parameters are the same of Fig. fy 

Notwithstanding the considered cellular flow has only 
one spatial scale, it is interesting to compare the results 
with the relation pit] 



Vf 



exp 



Un 



Vf 



(14) 



originally proposed by Yakhot ||35[ and Shivanshinsky 
p6[ for (multi-scale) turbulent flows, where U rms is the 
turbulent intensity (i.e., the root mean square velocity) 
and a = 2 and d = 1 are two parameters depending on 
the flow. Indeed, Eq. (|l4|) has been frequently used in lit- 
erature also for non turbulent flows and various values of 
a have been reported RpJ]- In particular, for the same 
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FIG. 4. The measured Vf/vo as a function of U/vp (□), 
data obtained using the measured values in Eq. (hj) with 
a — 2 and d — 1/2 (o) (see text), the function if)(U/vo) ( |ll| ) 
for f3 = 0.875 (solid line), and the upper and lower bounds, 
f3 — 1 in Eq. ( |ll| ) and Eq. (|l3|), respectively (dotted lines). 



A. Effect of small scales 



We consider now the following generalization of 



Eq. 



. /2tt \ /2tt 
u x (x, y) = U sin [ —x I cos I — y ) + 



N 



^ Qn sin \^-k n x + cos yj^ k ny 



ii , (.c. //) = - r cos ( ^-x j sin y-^y 



(15) 



N 



^Q„cos ( ^k n x + (j>* J sin (-j^ k ny 



where N is the number of different scales present in the 
flow, {fc„} are integers giving the ratio between the dif- 
ferent spatial scales, Q n is the velocity intensity at scale 
~ l/fc„, {4>n, <j>n} are (time-independent) phase differ- 
ences. 

In Fig. H we present two snapshots of the front for dif- 
ferent parameters values. By comparing with Fig. || it is 
clear the presence of small structures in the front due to 
smaller scales in u. 

We computed vt for N = 1 (two-scales flow) and N = 2 
(three-scales flow) with different values of k n , Q n and 



4 



— 1/3 

random phases. In the case N — 2, Q n — U k n ' 
has been chosen as a caricature of the power spectrum 
of three dimensional turbulence. The results, compared 
with the one-scale flow (Q n = 0) are summarized in Fig. 
where Vf is rescaled with vq and reported as a function of 
U Tms /v (U Tms = [7^(1 + Q 2 /[/ 2 )/2, with Q 2 = J2 n Ql)- 
As one can see, they roughly collapse onto a single curve 
together with the one scale results. 
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FIG. 5. Images of the burnt area (black) for U 
(above) and U = 10, Q = 2 (below) with k = 
vq — 0.5. The different scales are clearly visible. 



= 3 
where 



Let us conclude this section with some remarks. For 
the one scale flow in the steady case, we found that 
Vf ~ U/hiU as the flow intensity U (or U lms equiva- 
lently) is very large. This corresponds to the large U 
limit of the Yakhot like formula ( |l4| ) with a = 1. Our 
results indicate that the main effect of the small scales is 
to renormalize the stirring intensity. Indeed once taken 
into account U Ims all the curves roughly collapse (Fig. ||). 

However, because of the limited range of spatial scales 
here investigated it is difficult to say something definitive 
on the front propagation in a multi-scale velocity field. 

From Fig.s | and| one can sec that the introduction of 
small scales on u causes a roughening of the front shape 
at the same scales of the velocity field, and has just a mi- 
nor effect on the front speed. The problem of the effect of 
small scales on the front wrinkling has been studied also 
by Denet 6^] in the context of the flame propagation 
equation (FPE) (which is similar to the Kardar Parisi 
Zhang equation Pq|). In that model the small scales in- 
fluence on the large scale motion seems to be important. 
However, it is not possible to make a direct comparison 
between the two results, because the FPE and Eq. (||) 
generate different dynamics. 



IV. NON-STATIONARY CELLULAR FLOW 



The physical information one can extract from the 
above result is that the front speed is essentially deter- 
mined by the large scale behavior of the velocity field (i.e., 
the flow intensity). Indeed, as previously observed 
the absence of open channels can be more important than 
the detailed multi-scale properties of the flow. 
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FIG. 6. Front speeds Vf/vo as a function of ?7rms/i>o for 
the flow ( |l5| ) for various k n , Q n -values and vo values. Filled 
line is the one scale data with vo = 1.0. Two scales data, 
with v = 0.5, are: Qi = U/2 and ki = 3(A), Qi = U/5 
and fci = 3(x), Qi = U/2 and k% = 5(»), Qi = U/5 and 
ki = 5 (o). For a turbulent caricature field we used no = 0.5, 
Q„ — Ukn 1 ^ , k n = 3™. The reported data refer to 2 scales 
(□) and 3 scales (O), only. 



We now consider the problem of front propagation in 
the time dependent cellular flow 



u x (x,y,t) 



[/sin [f-x- 



Bsm(wt)] cos [^-y] 



,(x,y,t) = -[/cos [^x + B sin(wt)] sin [2Zy] 



(16) 



where the term B sin(w£) mimics lateral oscillations of 
the roll pattern, which are produced by the oscillatory 
instability |}0| . As one can see the steady case (0) corre- 
sponds to B = uj = Q. When B,oj ^ chaotic Lagrangian 
trajectories appear pO| , pO[ . The presence of complex par- 
ticle trajectories constitutes a step toward more realistic 
flows. 

We are mainly interested in addressing the two fol- 
lowing issues. First, since trajectories starting near the 
roll separatrices typically have a positive Lyapunov expo- 
nent, it is natural to wonder about the role of Lagrangian 
chaos on front propagation. Second, from previous works 
Ipgyioll , we know that for the time dependent flow (|l6| ) 
the transport properties are strongly enhanced, therefore 
it is worth to see if similar effects are reflected also in the 
front speed. 



A. Effects of chaos: transient dynamics 

A direct consequence of Lagrangian chaos is the expo- 
nential growth of passive scalar gradients and material 
lines 1 20 2^] : a (passive) material line of initial length £q 
for large times grows as 



5 



i(t) ~ £ e At . (17) 
Where A is the first generalized Lyapunov exponent, 

A= lim Inn WMV 

i—oo |5r(0)|-»0 i \ \Sr(Q)\ I 

which is in general larger than the maximum Lyapunov 
exponent [MpTI. 





FIG. 7. Snapshots at two successive times, t = 3.6 and 7.5, 
of the evolution of passive (top) and reactive line of material 
for two values of Vo (middle vq — 0.7 and bottom vo = 2.1) 
for U = 1.9, B — 1.1 and u) — 1.1U. The initial condition is 
a straight vertical line. 

The average in the previous equation is taken along 
the Lagrangian trajectories. In the presence of molecu- 
lar difTusivity, the exponential growth of l{t) stops due 
to diffusion |2(J and chaos is just a transient pi] . 



(due to stretching and folding), in the reactive one after 
a number of folding events structures on smaller scales 
are inhibited as a consequence of the Huygens dynam- 
ics: the interface between the two phases merges. This 
phe nomenon is responsible for the formation of pockets 
p5| , ^6 31 1 . Of course, "merging" is more and more effi- 
cient as vo increases (compare the middle and lower pic- 
tures of Fig. 0). 

In Fig. U we show the time evolution of the line length, 
C(t), as a function of t for the passive and reactive ma- 
terial at different values of vq. It is clear from the figure 
that while at small times both the passive and reactive 
scalar lines grow exponentially with a rate close to A, 
at large time t > t* (where t* is a transient time de- 
pending on vq) the reacting ones stops due to merging. 
At stationarity, the front length varies periodically with 
an average value depending on vq. A rough argument 
to estimate t*, is the following: two initially separated 
part of the line (e.g., originally at distance £q) become 
closer and closer, roughly as ~ £oexp(—At). When such 
a distance becomes of the order of v$t the merging takes 
place, hence to leading order 



t oc — In 

A \v 



(18) 



In the asymptotic state (t > t*) both the spatial and 
temporal structures of the flow become periodic. 



B. Effects of chaos: asymptotic dynamics 
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FIG. 8. C(t)/C(0) as a function of time for U= 1.9, B = 1.1 
and u)= Lit/ for the passive (+) and reactive case: from top 
v = 0.3 (x), 0.5 (*), 0.7 (□). The straight line indicates 
the curve exp(Ai) with A « 0.5, which has been directly mea- 
sured. 

For reacting scalars something very similar happens. 
Let us compare the evolution of material lines in the 
passive and reactive cases (see Fig. |?]). While in the pas- 
sive case structures on smaller and smaller scales develop 



Let us now switch to the effects of Lagrangian chaos on 
the asymptotic dynamics of front propagation. An im- 
mediate consequence of Eq. ( [18] ) is that the asymptotic 
front length (^) behaves as Lf ~ v^ 1 for enough small 
values of vq. Indeed, 



Lf ~ Le 



A/* 



L 2 A 



(19) 



which is in fairly good agreement with the simulations 
(see Fig. ||). It is interesting to note that in the steady 
case a different scaling can be seen. For the sake of com- 
pleteness, in the inset of Fig. 9 we display the dependence 
of the front speed on vq for both the time-dependent and 
time- independent flow. For very small vq, when chaos is 
effective in enhancing the front length, one observe an 
increasing of the front speed. At large values of vo the 
time-independent flow has a larger front speed than the 
time-dependent one. As we will see in the next subsec- 
tion, this is a consequence of the mode-locking of the 
dynamics, which maintains constant the value of Vf. 

It is worth remarking that even if the scaling jl9| ) 
holds when chaos is present, in general it is not pecu- 
liar of chaotic flows. For instance, for the shear flow 
(u x = U sin(y) , u y — 0) one has Vt = U + vq. On the 
other hand from Eq. (0) Vf ~ L(Vq. Therefore, even if the 



G 



shear flow is not chaotic Lf ~ 1/vq for U/vq ^> 1. From 
the previous discussion, it seems that the front length 
dependence on vq is not an unambiguous effect of chaos 
on the asymptotic dynamics. But, comparing Fig. 2 with 
Fig. 7 appears that the spatial "complexity" of the front 
in the presence of Lagrangian chaos is higher than in its 
absence. 
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FIG. 9. The average front length Lf as a function of v 
for the time dependent flow (o), with U — 1.9, B = 1.1, 
uj = 1.1U and the time independent case (x) with U = 1.9. 
The straight line indicates the 1/vo behavior. In the inset it 
is displayed Vf versus Vo for the time dependent (□) and for 
the time independent case. 




FIG. 10. Lf/Wf as a function of no for the time dependent 
(o) and independent (x) cases for the same parameters of 
Fig. 9. 

In the sequel we introduce an indicator to quantita- 
tively evaluate this qualitative observation. 

Let us call Wf the size of the region in which burnt 
and unburnt material coexist. Introducing a measure, 
in that region we can defined Wf as the standard 



deviation of n(x) |24 



\dj(x)\ 
J dx\dj(x)\ 



(20) 



where 6{x) = l/L J Q L 9{x, y)dy, i.e. 



W f 



x /i(x)dx — ( / x/i(x)dx) 



1/2 



(21) 



For a simple shear flow Wf and Lf display the same kind 
of dependence on vo (actually they are proportional). In 
generic chaotic flows there is an increasing of the front 
length, while chaotic mixing induces a decrease of Wf. 
This is indeed what one observes in Fig. 10, where we 
show the ratio Lf/Wf both for the non-chaotic and the 
chaotic flow. For the latter this ratio diverges for very 
small vo values as a signature of chaos. From a physical 
point of view the ratio Lf/Wf is an indicator of the spatial 
complexity of the front. Indeed it indicates the degree of 
wrinkling of the front with respect to the size of the re- 
gion in which the front is present. Loosely speaking, we 
can say that the temporal complexity of Lagrangian tra- 
jectories converts in the spatial complexity of the front. 



C. Front speed dependence on the frequency 

For passive particles transport in the flow (^), it has 
been found that the eddy diffusivity coefficient D c ff(u) 
displays a complex behavior with resonant-like patterns 
[^9|,^0|, in which the values of D e ff(u>) are orders of mag- 
nitude larger than the stationary-flow value, -D c ff(0). The 
physical mechanism responsible for the resonances is re- 
lated to the interplay between the oscillation of the sep- 
aratrices and the circulation inside the cell. When circu- 
lation and oscillation "synchronize" , a very efficient and 
coherent way of transferring particles from one cell to 
the other takes place. Does it happen something simi- 
lar to the front speed in the reactive case? In Fig. |ll| 
we report Wf(e) as a function of e = lo/U. As one can 
see, Vf(e) varies both above and below the time inde- 
pendent value, Vf(0), and its range of variability (about 
30% above or below the time independent value Wf(0)) 
is very small compared with that of the diffusion coeffi- 
cient. This can be understood as a consequence of the 
inequality Vf < U+vq. If one changes the value of vo then 
the curve display a small shift up (increasing vo) or down 
(decreasing vq) with some slight variations of the values 
of the peaks. This behavior is very robust with respect 
to changes of the flow parameters (U and B). Moreover, 
Vf (e) is a piece-wise linear function of e with slope given 
by U multiplied by a rational number. This phenomenon 
-called frequency locking ^,^3|- naturally arises in some 
non linear oscillators which are periodically forced, when 
the oscillator synchronizes its frequency to the one of the 
periodic forcing. 
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The frequency locking of the front speed can be un- 
derstood with the following argument. At large times, 
t > t* , the front is time and space periodic. This means 
that after a time period T, the front is rigidly trans- 
lated in the ir-direction by S, which is the spatial period. 
Due to the spatial periodicity of the flow ( Jig ) S = 2nN 
(where TV integer). Now, as confirmed by the simula- 
tions, one expects T to be a multiple of the oscillation 
period T a = 2tt/uj so that: T = MT a (with M integer). 
On the other hand, the front speed is nothing but S/T 
so that 



not observed because of the enhanced transport proper- 
ties of the flow Fig. |ll|). 



V{ 



S 2ttN 



N N 
— u> = — Ue , 

M M ' 



T MT 

which is indeed the behavior we observed in Fig. Ell 



(22) 



1.2 
1.1 
1 

S 0.9 
> 

0.8 
0.7 
0.6 



1(2,1) (1.1) 



(1,2) 



(1,3) 



(1,4) 



0.5 



1.5 



(1,5) 



' (1,6) 



2.5 



FIG. 11. Vf(e) as a function of e = lo/U, for the flow ( pi| ) 
with U — 1.9, vq — 0.2 and B — 1.1. The straight horizontal 
line indicates the front speed for the steady case, B = 0. The 
dashed lines indicates the curves U eN/M for different (N, M) 
integers. 

Varying uj the periods S = 2ixM and T = NT a with 
a given M and N can loose their stability so that a new 
couple of N, M values is selected. This explains the pres- 
ence of different linear behaviors. By generalizing the 
one dimensional model ( |l0| ) to the time dependent case 
one can qualitatively reproduce the behavior of the front 
speed dependence on the frequency (see Appendix |b|) . 

It is interesting to compare the behavior of front propa- 
gation with the time dependence (|l6| ) with previous stud- 
ies that considered a different time dependence, i.e., with 
U -> Ucos(ujt) and B = in (|l|). With this choice the 
overall effect (as recognized by most of authors [ ^T|j42| 
see also [[h|) is a depiction of the front wrinkling at in- 
creasing the flow frequency. As a consequence, a strong 
bending of the front speed with respect to the steady case 
has been observed. This phenomenon has been quanti- 
tatively understood for the case of time dependent shear 
flows (i.e., u = ([/sin(wi) sin(y), 0]) by Majda and col- 
laborators . With the choice (|lq) such a depletion is 



V. FINAL REMARKS 



In this paper we studied thin front propagation in 
steady and unsteady cellular flows. In particular, we in- 
vestigated the behavior of the front speed and the front 
spatial structure at varying the system parameters. 

As far as the one-scale steady case is concerned, we 
were able to give a quantitative description of the front 
speed by means of a simple one dimensional model. For 
large flow intensity U (or U rms equivalently), the front 
speed behaves as V{ ~ U ims / \og(U rms ), which corre- 
sponds to the asymptotic behavior of the Yakhot-like for- 
mula ([li]) with a — 1. Moreover, small scales structures 
have been added to the flow in order to study their effect 
on the front speed. Numerical simulations show that, 
once Vf is rescaled with vq and reported as a function of 
U lms /vo, the results for the one-scale flow and those for 
more than one-scale fairly collapse onto a single curve. 
Therefore, the front speed is essentially determined by 
the large scale behavior of the velocity field. 

Small scales spatial structures may also be induced by 
Lagrangian chaos. In this respect, on the basis of our 
results on the unsteady cellular flow it is interesting to 
remark that the effect of chaos is limited to a transient, 
in which the front behavior is close to the passive scalar 
case. Asymptotically, the reacting term induces a dras- 
tic regularization on the front evolution, suppressing the 
effects of small scales and Lagrangian chaos. Indeed the 
front propagates periodically and displays a frequency 
locking phenomenon. 

The only asymptotic effect of Lagrangian chaos is in 
the structure of the front which is more and more wrin- 
kled as vo approaches to zero. On the contrary in the case 
of steady velocity fields (regular Lagrangian motion) the 
degree of irregularity does not change with vq. As an in- 
dicator of the spatial "complexity" of the front we used 
the ratio between the front length and the width, Lf/Wf 
which is large (diverging as vq — * 0) for the unsteady case 
and is roughly constant for the steady one. 
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APPENDIX A: NUMERICAL ALGORITHM 



APPENDIX B: FRONT SPEED LOCKING 



In numerical approaches one is forced to discretize both 
space and time. We introduce a lattice of mesh size 
Ax and Ay (where we assume Ax — Ay) so that the 
scalar field is defined on the points x„ j7n = (nAx,mAy): 
6 n ,m(t) = 9(nAx,mAy,t). 

The time discretization implies a discretization of the dy- 
namics. Looking at the G-equation (|3|) one immediately 
recognize two different terms: the advection term u-VG, 
accounting for the transport properties of the flow, and 
the "optical" term vo|VG|, which locally propagates the 
front in a direction perpendicular to it with a bare veloc- 
ity v . 

Let us call F At the Lagrangian propagator for the dis- 
cretized advection equation, <9jx = u • VG. Then, given 
the field at time t one can computes the field at time 
t + At using a two steps algorithm: 



using the Lagrangian propagator, F At (x), one 
evolves each point of the interface between burnt 
and unburnt region; 



2) at each point of the evolved interface one constructs 
a circle of radius vq At, obtaining the new frontier 
as the envelope of the circles. 

t __ Af t+At 
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FIG. 12. Pictorial scheme of the numerical algorithm for 
the geometrical optics limit. 

The numerical algorithm can be easily implemented us- 
ing a time reversal procedure: starting on a grid point, 
x n ,m, of the scalar field at time t + At one applies the 
backward evolution obtaining the point y = F _At x n)m 
at time t. Around y we construct the circle of radius 
vo At. If in this circle there is at least one burnt point 
of the scalar field at time t, we fix 9(x n>m ;t + At) = 1 
otherwise 0(x n>m ; t + At) = 0. 

It is worth to note that one has to care about the radius 
of circle v$ At that has to be larger than the grid-size Ax 
(empirically we found that v^At > 3 — 4 Ax it is enough 
to give good results). 

Typical values used in our simulation are At — 0.02, 
Ax = 27r/2048. The backward Lagrangian integration 
F~ A *(x) has been performed with a 4-th order Runge- 
Kutta algorithm. 



Frequency Locking arises in many physical systems 
ranging from Josephson-junction arrays to chemical reac- 
tions and non linear oscillators p2]-p4]] . The basic mecha- 
nism is a resonance effect between two oscillators or when 
an oscillator is coupled with an external periodic forcing. 
In the last case the system synchronizes with the exter- 
nal forcing making its internal frequency commensurable 
with the external one. Almost all the systems displaying 
frequency locking can be mapped to the damped forced 
non linear oscillators [E2|: 



#0 



rlf) 

+ j3—+~fsm(0)=5 + acos(wt). (Bl) 
at 



The solution 9(t) is periodic and the frequency, i.e., the 
average angular velocity, turns out to be: 



y dt> 



9 M 
= hm - = —uj 

t^oo t N 



(B2) 



with M,N integers. Moreover, if (B2) is realized for a 
certain set of the parameters there always exist an entire 
interval around their values where (B2) holds with the 



same values of M and N. This kind of behavior persists 
also when a = and for other kind of non linear terms 
(i.e., the third term of the l.h.s.). An exhaustive descrip- 
tion of such a phenomenon can be found in Refs. p3|. 




FIG. 13. vm(c) ad a function of e = u/U, for the model 



(B3) with U = 1.9, v = 0.2 and B = 1.1. The dashed lines 



indicates the curves UeN/M for different N,M integers. 

Coming back to our system, we can generalize the ld- 
model to the time dependent case: 



dx M (t) 
dt 



= v + U\ sxn(xM + B sin(wt))| . 



(B3) 



Note that in principle one should also take into account 
the dynamics of yM, but for the sake of simplicity we 







present just the y-independent version of the model. This 
model, although very idealized, is able to reproduce be- 
haviors qualitatively similar to the ones observed in the 
simulations (compare Fig. |l3|) with Fig. B y changing 
of variables z(t) = Xm(P) +Bsin(uit) Eq. (BS) reduces to 



dz(t) 
dt 



— v + U\ sin(z))| + Btvcos(u)t) . 



(B4) 



which corresponds to (Bl) for a — 0, and for which fre- 
quency locking has been studied in details. It is inter- 
esting to quote a recent work |Q which has studied the 



problem of locking in a model very similar to (B4) but 
in presence of noise. They found that the locking phe- 
nomenon is rather robust under the effect of noise and, 
moreover, it gives rise to resonances in the diffusion co- 
efficient. All these results are qualitatively very similar 
to the behavior of the system here studied. 
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